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Abstract: Consider the classical XY model in a weak random external field pointing along 
the Y axis with strength e. We prove that the model defined on with nearest neighbor 
coupling exhibits residual magnetic order in the horizontal direction for arbitrarily weak 
random field strengths and, depending on field strength, sufficiently low temperature. 



1. Introduction 

In this paper we study an interesting phenomena which leads to ordering at low tem- 
peratures for spin systems with continuous symmetries: random field induced ordering 
(RFC). To fix ideas and introduce the model which is studied below, consider a classical 
0(2) model on Z**. For each x e X'^ let ax E S^. A vector of spins a = (c^x)xezrf will 
be called a spin configuration. Let {'>^x{(^))xgZ'I auxiliary i.i.d. family of standard 

normal random variables, with co representing an element of an auxiliary sample space 
n on which the Kx's are defined. For any fixed spin configuration and any bounded 
region A C Z'*, we can then define the (random) Hamiltonian via 

-n%{a\iT') = -^ X] [crx-(ry]^ + ej^ax{cv)e2-crx (1.1) 

(x-i/)nAjt0 xeA 

where we set ax = ax for x G A"^ and allow ax G to be arbitrary within A. Here (xy) 
indicates that x, y are nearest neighbors with respect to the usual graph structure on Z'^ 
and e-[,e2 denote the standard orhonormal basis of K^. Defining finite volume Gibbs 
measures by 

and denoting the corresponding Gibbs state by (•) = (■)^, the question is whether and 
when, in terms of j6, e and ctq, residual magnetic ordering occurs in the limit as A f 
T.^. Below we discuss in more detail previous work on the subject, but prior to the 
present work it was expected that ordering does occur in dimension > 3 while in 
dimension d = 2 there was uncertainty about the low temperature behavior [1, 8, 9, 29]. 
In general, when ordering does occur, the ordering is expected in the horizontal direction 
±ei whenever the projection a^ ■ e^ is either uniformly positive or uniformly negative 
on Z'*. In this paper we demonstrate that ordering occurs in the ±ei direction only if 
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d = 3 and only if the boundary condition = ±ei. Our framework can be extended to 
prove the same statement for the two dimensional system and also to treat more general 
boundary conditions = u for some u G S^. Because these statements increase the 
technical complexity substantially without shedding further light on the methods, we 
leave them for future work. When possible, we will comment on modifications which 
are needed to treat these extensions. 

Despite the fact that we only provide a proof in three dimensions, the result is most 
interesting if = 2. To explain why, we recall the behavior of related models. 

• The Pure 0(2) Model. In this case we take the same setup as (1.1) except that we 
set e = 0. When d = 2 there is no residual magnetic order in the thermodynamic 
limit (this is the content of the Mermin- Wagner theorem, see [23, 22, 7, 16] among 
many other works). There is however a Kosterlitz-Thouless phase transition [17] 
expressed by a change in the behavior of the decay of the spin-spin correlation 
function {dx ■ dy) with \x — y\. If d > 3, residual magnetic ordering occurs [18]. 

• The Random Field Ising Model. In this case we constrain the spins to be ±62/ 
replace di/ by unbiased counting measure and otherwise retain the setup of (1.1). 
When d > 3, residual magnetic ordering occurs for ±62 boundary conditions. 
This was proved for ground states in a system with weak disorder in [20] and 
for low temperature and weak disorder in [4]. On the other hand in dimension 
d = 2 it was proved in [2], see also [21], that there is a unique infinite volume 
Gibbs state at arbitrary strength e of the disorder. 

• The Random Field Gaussian Model. In this example we replace the vector val- 
ued spins (Tx S by a scalar field (px G K, replacing 62 by 1, and replacing dv 
by Lebesgue measure on ]R. The model appears in [10], see also [27] for related 
work. Because the underlying configuration space is no longer compact, the ex- 
istence of infinite volume Gibbs states is a (somewhat) nontrivial issue and, to an 
extent, is the analog of the question of residual ordering for compact spin spaces 
with a continuous symmetry. 

To get an indication as to what should be expected we compute two point cor- 
relations in finite volume. Fixing a finite volume A, let — A^^^ denote the discrete 
Laplace operator on A with Dirichlet boundary conditions. Taking boundary 
conditions (px = and averaging over the ocz's gives 

n{Mj)t"'] = -A;,i(x,y) +e2A-2(x,y). 

If e is 0, these correlations are finite (uniformly in A) only if d > 3. In this case one 
can define infinite volume Gibbs states for the field (px- When d = 2, one must 
instead view the infinite volume measure as existing on the space of gradients. If 
6 > and if d < A, the second term on the RHS grows with A while, when d = 2, 
it grows even after taking gradients in both arguments x, y. 

• Order-by-Disorder. On the other hand, the phenomenon of "Order-by-Disorder" 
provides related examples of systems which exhibit ordering due to various 
types of fluctuations. The most relevant example, first considered by Henley 
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[19], concerns a model Hamiltonian on 1?- of the form 

\V-v\V=^ \\^-v\\ 

with < 2/2. The ground-states for this (frustrated) system are obtained by 
choosing a purely anti-aligned configuration of spins on each of the even and odd 
sub-lattices of and are thus parameterized by two angles: an angle between 
the spin at (0, 0) and the ei-axis and relative angle between the spin at (0, 1) and 
(0,0). The degeneracy of ground-states is partially lifted under the introduction 
of two types of "disorder". The first type consists in passing from to positive 
temperature, see [3] for mathematical justification of the effect in this case. More 
relevant to the RFIO is a second mechanism: site dilution. Vertices of 1?- are 
deleted from the system independently with probability p ^ 1. According to 
the calculations in [19], at and low temperature the system prefers the ground 
states with relative angle between spins at (0, 0) and (0, 1) to be fixed at ±^/2. 

Besides the obvious differences between this setup and ours, there is one cru- 
cial technical difference. In the site diluted model, large fluctuations due to ran- 
domness are substantially weaker than those encountered in the analysis of the 
RFO(2) model. In particular, there is an analog to the field gx introduced below 
but the fluctuations of this field are about as singular as the four dimensional 
version of the RFO(2) model. While the site diluted order-by-disorder problem 
has not been rigorously addressed, this feature suggests the conclusions in [19] 
are reliable. It would be interesting to see if our methods can be adapted to this 
case. 

Of the first three examples, the most worrying from the perspective of proving mag- 
netic ordering in the RFO(2) model is the last one: If we represent spins via angular 
variables - = (cos(0v)/ sin(0;i)) - and make an expansion of (1.1) around Q-^ = (or 
any fixed angle i/^ for that matter) our model looks exactly like the random field Gauss- 
ian model. This suggests that the whole ansatz of having order may be flawed since 
fluctuations in the latter model are so strong. 

Beginning in the early 1980's, RFIO was the subject of a number of publications in 
the physics literature, see [8, 9, 24]. One group [8] concluded there is a low tempera- 
ture paramagnetic phase The other [24] concluded there is an intermediate-temperature 
ordered phase from which they extrapolate the low-temperature behavior. Interesting 
tangentially related rigorous work was done in the 1990's on ground states in the strong 
field regime in [14, 15]. 

The question of ordering in the RFO(2) model has also appeared in recent literature. 
For Bose-Einstein condensates in optical traps, the effect was suggested as a possible 
response to the presence of certain kinds of experimentally realizable disorder [29, 26]. 
Here the (pseudo-)spin variables arise from internal structure of the atoms, the tuning of 
interactions and the structure of the optical lattice. This type of phenomenon appeared 
as a possible mechanism for the splitting of Landau level degeneracy in graphene [1]. Fi- 
nally, van Enter and coauthors [11, 12, 13] came to this question during investigations of 
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whether the spatial Markov property which characterizes Gibbs measures is preserved 
under various coarse-graining procedures. 

Little mathematically rigorous progress has occurred for RFIO except regarding qual- 
itative ground state behavior and the mean field approximation [11, 12, 13, 28, 29]. Most 
recently, the author derived results [5] consistent with the picture presented above for 
the RFO(2) model with a Kac potential. The fundamental limitation of that work is that 
the results are only valid if the range of the interaction potential is taken to diverge as 
a polynomial in e^^. This limitation is not entirely technical as we show below. Even 
the nearest neighbor RFO(2) model has a fundamental length scale which is set by the 
behavior of energetics. 

1.1 Main Result. The main theorem we aim to prove is the following: Let e be fixed, 

L ~ e^^ log^£ 

so that L = 2^ for some k G N, and 

Qi{z)=z + {d,l...,L-lY. 

A subset A of "ZJ^ will be said to be L-measurable if A is a union of blocks QL{r) so that 
r G LZ'^. We define block average magnetizations 

Theorem 1.1 Let d = 3. For any ^ G (0,1) small enough and any e G (0,eo(^)) there is 
/3o(e, ^) so that if fi > fio the following holds: Suppose that Aj^ is a sequence of L-measurable 
volumes which increase to 1? in a van Hove sense. For almost every a; G Q, there exists an 
L-measurable subset T)co C Z*^ and an Nq{lo) G N such that: 

(1) iD^^nAwl < Ce-^|i°s^l'|AN| ifd = 3,forallN>No{w). 

(2) For each z G Ajv with Ql{z) n D^; = and 

ii(Mz)r-^iib<?- 

Here C, c are universal constants. 

Looking into the proof, we have a bound on the transition temperature of j6(e) < 
e^'^Poly{ \ loge|) (we get this sort of bound also in the two dimensional case). This is 
inline with the predictions of mean field theory, see for example [?]. On the other hand, 
the analysis in [24] suggests a random field strength independent transition temperature 
in the two dimensional case. Whether this prediction is correct and what might take its 
place if it is not correct is unclear at present. Extensions. We have already mentioned 
that Theorem 1.1 may be extended to include d = 2 and also other (constant) boundary 
conditions. We emphasize in the latter case that the derived conclusions will be the same 
as in Theorem 1.1: if the constant boundary condition u satisfies u ■ e^ > 0, then deep in- 
side the volume Mz will still be very close to ei . Another direction in which the theorem 
can be extended is to incorporate spins cr^- G S". In this case we take the disorder a^to 
be n — 1 dimensional i.i.d. standard Gaussian vectors supported on some hyperplane 
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through the origin. Ordering then occurs perpendicular to this hyperplane. Finally, 
Gaussian disorder is taken only for convenience and Theorem 1.1 can be extended to 
incorporate any choice of disorder for which the estimates of Lemma 11.1 (l)-(4) hold 
(Statement (5) of that lemma follows from Statements (1) and (3)). 



1.2 Some Intuition and an Outline of the Paper. In this subsection we give the heuris- 
tics and basic method of proof which lead to our result. The discussion will be given 
here for both d G {2, 3}. Here is the Kronecker it is 1 if d = 2 and otherwise. This 
explanation is a bit impressionistic to avoid some technical details. 

First, let us indicate why one should expect ordering to occur in the horizontal direc- 
tion (if there is to be ordering at all). We restrict attention to boxes C of sidelength 

i < € \ log e I ~ for numerous reasons which we hope will become apparent as the 
discussion proceeds. 

Let Cx be given in polar coordinates by an angle 9x . We perform a "spin wave" analysis 
supposing that in Qf 6x = ip + 9x with 9x small and xp fixed. Then we expand —'Hq,- (with 
free boundary conditions) in dx variables. Keeping only terms up to second order in dx 
and e we find 

^2 

sup --Hq,(0) = --cos2(i/^) dcx^-^ ■Oix + 0{e\ 

^ V ■' 

where A is the Neumann Laplacian for Q( and dcx = oCx — \Qe\^^ ^zeQt direc- 
tions of presumed ordering are obtained by optimizing the first term in xp and ignoring 
Term L An important additional observation is that the optimal choice for the deviation 
variables 6x is 6x = cos{xp)g^ where 

:= -eA-'&x. 

In evaluating the validity of this computation we run into some constraints on the box 
size £. Typically, 

Y cixA-^ ■ ocx ~ e^e^ log-^^'"* i. (1.2) 

So, for d > 3, the central limit theorem implies that Term I will be of lower order 

2 

(with high probability) if £ 3> e^'^. For d = 2 there is crucially an extra factor of 
log£ in (1.2). This allows us, for appropriate choice of 7 < V2, to suppress Term I for 
i ^ e^^\ log e I ^ 5 ^ with probability exponentially small in — | log e | ^ ^^"^ . Together these 
considerations give a lower bound on £ for d > 2. 

On the other hand, in low dimension the behavior of the field provides an upper 
bound on £. This field has a typical order of magnitude in dimensions d = 2,3,4, oi 
i, y^logT respectively. Thus the calculation cannot be taken too seriously for large 
boxes as in this case the maximizer is inconsistent with starting assumption that Qx is 
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(uniformly) small. In two and three dimensions, we arrive at the following constraints: 

e"^|loge|"2"T'<£<e"^|loge|"5 ifrf = 2, (1.3) 
< £ < e-2|ioge|-2 ifd = 3. (1.4) 

where the | loge|^2 factor on the right hand bound was inserted to account for possible 
fluctuations of | j 1 1 oo . This indicates the marginal nature of the two dimensional model. 

From a technical perspective, the spin wave calculation seems to depend on uniform 
control of 6x- Unfortunately, the only control we have, stated in Proposition 8.1, concerns 
the Dirichlet energy of spin configurations in Qf, <fQ( (cr) := L(j,,/)cQ(['^t — c^/]^- It is 
energetically favorable for 

< epilog el-^^-^^rf (1.5) 

but no better. Indeed, from the spin wave heuristic our best guess on the maximizer of 
— "Hq^ is which has Dirichlet energy on this order. 

Because we can only control the oscillations of a in H^, we needed a different com- 
putational device which gives the same conclusions as the spin wave analysis. This is 
presented in detail in Lemma 8.2. The physical idea which sits in the background is that 
low energy spin configurations in the RFO(2) model consist of two superposed modes. 
In angular variables, with Cx replaced by the angle 6x, Ox consists of a slowly varying con- 
figuration (px, whose Dirichlet energy is negligibly small, and a quickly varying mode 
which is proportional to . A convenient way to distinguish between these two modes 
is to define (px by 

(Px = ex-cos{ex)g^. (1.6) 

If we regard this as a change of variable 9x ^ (px the Hamiltonian —"Hq^ transforms into 
a new energy functional 

^QfX'Px) = E -'Py)-'^ +jT,^nxCos^{(px) 

^ V ' 

II 

with 

y:\x-y\=l 

There are errors made in this transformation but we can control these errors using Sq^ {a) 
and the behavior of the two fields g^ and (Xx in Qg so long as £ is neither too small or 
too big as indicated above. Note that the potential ~ logej''^." typically. Thus if 
(px only varies at a much longer length scale than g^ we may pretend that has been 
replaced by this constant. Because Term II punishes configurations for which (px varies 
quickly, the conclusions of the spin wave computation are still valid. 

Using the above analysis, we next setup a coarse-grained Peierls argument using the 
scale £ as the microscopic scale and show how to extract energetic cost from the oc- 
currence of contours. The technical problem in implementing this idea is that the spin 
space is continuous and there is no microscopic surface tension. This is analogous to the 
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problem encountered in constructing a Peierls estimate for a "soft" Ising model in which 
continuous scalar spins sit in a double well potential. Limitations inherent in the bound 
(1.5) make the present analysis more challenging. To solve this problem, we take inspi- 
ration from previous work on Kac models, see [25] and also [5]. This method requires 
us to use two scales; the scale £ and a second scale L ^ £ with contours defined relative 
to the second scale L 

Before describing contours let us fix the scales £,L. A bit later we will see that it is 
important that £ < £''^\loge\'^^'''^^. We take 

£ ~ e^^l loge|^^ and L ~ loge|* if = 3, if d = 3, 

£ £-^\loge\-''^-'^^^^ and L ~ e^^l loge|-''2+'^™ if = 2. 

To define contours, we first classify cubes Qf as either good or bad relative to a spin 
configuration a. Fixing cr, call a box Q{ of side-length £ bad for a if either the Dirichlet 
energy £Q^{cr) is substantially larger than 4e'^£'* log''^'' £ or if the average of spins in Qi, 
o'iQi) '■= ^ ^xGQt ^x, is far from ±ei; see § 3 for further details. A semi-precise defi- 
nition of a contour F is as follows: contours are maximally connected clusters of boxes 
Qi of side length L so that within distance 2L of Qi there is a cube Q{ of side-length £ 
which is bad for cr. 

The reason one introduces the second scale L is that it allows us to do surgery at the 
boundary of a contour to artificially create surface tension. Given a spin configuration cr 
and an associated contour F let (5(r) denote the neighborhood of radius L around F. The 
idea is to compare a with a new spin configuration a which agrees with either cr or the 
reflection of cr across the £2 slxis on each component of A\S (F) . We construct a to have an 
angle uniformly close to either or zr on the whole of F. The skeleton of this argument 
is given in §§ 5, 6 and 7, proofs of lemmas from these sections appear in §§ 10.1, 10.2 
and 10.3. 

The key, and most labor intensive, part of the construction requires us to do surgery on 
the spin configuration cr at the boundary of the contour, i.e. on (5(F)\F. We will assume 
the randomness is "uniformly typical" in the discussion of this procedure. Allowances 
must be made in the actual analysis and precise conditions we require of the randomness 
on a contour F are detailed in § 3.1. Not all contours satisfy these conditions, but we state 
lemmas which quantify the size and sparsity of regions where our requirements fail. 
The probabilistic bounds required to prove these lemmas and the proofs of the lemmas 
themselves are postponed until § 11. 

By definition, cr is good in (5(F) \F: for each cube Qg within Lof^(F)\F 

£Q,{cT)<4e^£'\og'^'^£, 

and £^'^ HxgQi '^^ close to ±ei. Combining these to facts, one can show that the sign 
±1 is constant over connected components of (5(F) \F. 

Fix a component R of ^(F)\F and suppose £^'^YiiX(^Q;(^x is close to ei throughout R 
(the other case is similar). The reason our analysis is restricted to low dimension is that 
if (i G {2,3}, the a priori bounds on £Q^{a) allow us to find a small neighborhood 
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R' D K SO that for vertices x in the graph boimdary of R' , 9^ G {—"/gf ^/s)- This point is 
addressed in § 9.1. 

Next, if we start with such a boundary condition ao for R' , we can look for maximizers 
of — 'Hr/((7|(7o). Actually, it is better to use Lemma 8.2 and look for maximizers after a 
change of variables. The idea is the same as in (1.6) but we use the Dirichlet Laplacian 
A^, in R' instead of the Neumann Laplacian. Set 

= e[-A^,+A]-V 

for all X in K'. The mass A = L^'^log^ L if d = 3. It is a technical convenience used to 
keep small in absolute value throughout R' (this is a point where one has to make 
modifications in two dimensions). We make the change of variables 

0, ^ 0, - cos(0,)^^,, =: (px (1.7) 

Because of the Dirichlet boundary conditions the restriction of cr to R'^, cr\j^,c, is invariant 
under this transformation. On the other hand, —'HR'{cr\(7o) transforms into 

^R'(fi#R'0 = Li cos(fv - 4)y) - 1 + jY^^nx COS^{(px) 

{xy}nR'^0 X 

where 

y-\x-y\=i 

This change of variables potentially costs us energetically, but we will more than offset 
any loss by what we gain in comparing a ioa (which is still to be constructed). 

One can show that if the randomness is well enough behaved in R' , there is a unique 
maximizer of ATj^/ and that this maximizer also has (px G ( — 7r/8, 7r/8) throughout R' . 
Because of this and the fact that typically nix ~ logej'^''^, one should morally regard 
the maximizer as behaving like a solution to the discrete elliptic PDE with mass 

- ^(px + I \o^e\^^-'''(px = 0. (1.8) 

This leads us to the separation of length scales at e | log e| t~ which we chose above. 
A solution to (1.8) is very close to deep inside R' only if the inner radius of R' is much 

larger than e \ log e| 2 . This issue is addressed in § 9.3. On the other hand, in order 
to distinguish between the pure phases ±ei, we need £ to be smaller than this scale since 

1 ''2''' 

I log e\~ will be the typical thickness of interfaces between pure phases. 

After replacing a by the maximizer of fC^i inside R' for each component R of (5(r)\r, 
we can invert the change of variables (1.7) to produce a new spin configuration , which 
disagrees with a only on (5(r)\r. Because is very close to one of ±ei deep inside each 
thickened boundary component of T, at negligible further cost we can force it to be 
exactly ±ei in the middle of each of the boundary components, calling the result cr^. 

Next we modify into a new configuration on . We keep fixed on the 
exterior component S of A]v\r and note the value sgn(c72 • ei) on fi ^(r)\r. We make 
the following transformation on the interior components Xj. If on Xj n 5 (T) \r sgn((72 ■ ei ) 
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agrees with its value on the boundary of the exterior component, we keep cr2 fixed on 
the entire interior component. Otherwise, we reflect (72 across the 62 axis on all of Ij. We 
call the result, which is defined only in Aj^\T, cr^. 

Let S{T) be the neighborhood of radius V2 of Y. To obtain a from cP, we construct 
(almost-)maximizers for — 'Hj(p)(-) (here we use free boundary conditions). This gives 
two configurations (7^ uniformly close to ±ei throughout S{T). The configuration a 
is chosen to agree with on S{TY and to agree with the appropriate maximizer of 
~^<5(r)(') inside ST. After accounting for all errors and also what we gain in the com- 
parison between a and a we are able to show 

-nAM^i) + nAM^i) > ce2|ioge|-24|r| 

for some C > 0. With this result in hand we can construct a Peierls argument and prove 
Theorem 1.1. 

The remainder of this paper is organized as follows. In § 2 we collect basic notation 
used throughout the rest of the article. § 3 defines what are contours in our model in 
a careful way, while § 3.1 addresses what sort of requirements we need from the ran- 
domness ctx in order for a subregion of Z'^ to be well behaved. Regions satisfying these 
requirements are called clean regions and in Proposition 3.3 and Lemma 3.4 we quan- 
tify the likelihood and sparsity of dirty regions. Proofs related to § 3.1 are postponed to 
§ 11. In § 4 we state the Peierls estimate. Lemma 4.1, and prove Theorem 1.1 on its basis. 
Assuming that a given contour is clean, § 5 outlines the construction of almost maximiz- 
ers to — while § 6 outlines the construction of ai above. In § 7 we construct a 
and then prove the Peierls estimate. Proofs of major estimates required in §§ 5, 6 and 7 
are postponed until §§ 10.1, 10.2 and 10.3 respectively. § 8, while technical, is key to the 
whole paper as it contains Lemma 8.2. Finally § 9 contains various technical lemmas 
needed elsewhere, for example in the analysis of maximizers of /Cr/ (^x|^r'0- 



2. Preliminaries 

Throughout Cx G S^. We will use the notation 61,62 for the usual orthonormal basis 
of in which sits. We will also use {ei, . . . , e^} to denote the standard basis in 
and R'*, the " being used to distinguish unit vectors which live in the spatial lattice from 



those which live in spin space. In general 
and / : ^ we let 



denotes the £P norm in R'*. If R C Z"^ 



\p,R 



L 



1/p 



We write x 



y x,y E'ZI^ a \\x 
d'R = 
d°R = 



- y\\i = 1. For any finite subset R C Z**, let 
{xeR: disti(x,R^) < 1} 
{xeR' : disti(x,R) < 1} 
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with distp {A,B) the Hausdorff distance between to sets in £P . | R | will denote (depending 
on the context) the cardinality (resp. volume) of a finite set R C Z'' (resp. bounded 
domain in R*^). 

The proof requires that we work on multiple scales. This is conveniently implemented 
by considering L' G {2^)i^^f^. We will use the notation 

Ql'(^o) = {z G Z'^ : z = xo + u for some t7 G {0,. . .,L' - l^} 

to denote an arbitrary box of side-length L' indexed (in the "lower left corner") by some 
fixed xo G Z'*. Also, let Bl'(^o) = xq + {-L' /2, . . . ,L' /2 - lY denote the cube with 
side-length L' centered (roughly) at xq. We shorthand these by suppressing reference to 
Xq, writing Qii, Bii instead when no confusion should arise. 

Let ccx, X G Z** be a field of i.i.d. Gaussian random variables with mean and variance 
1. We treat in full only the case of Gaussian disorder for now. Our present proof should 
work for any random variables with subGaussian tails by extending the estimates of 
§ 11. To what extent we need subGaussian tails, and not just exponential moments in a 
neighborhood of is not clear. 

We let Sa = [S^]^ for A C Z'^ denote the configuration space for a spin system 
restricted to A, with S = S^d. If Ai C A2 and a G then (j\ai G Sai denotes the 
natural restriction to (T" to Ai. Let 

(!■!/) 
x,yGR 

Given (Tq ^ S and a ^ Sr we may extend a to an element of S via cr\Rc = (JqIr^- In this 
case we let 

1 1 

- nR{o-\ao) = --£R{a) - - ^ {a^c- (^yf + ^ ^ OL^ei ■ (Tx- (2.1) 

We sometimes consider regions with free boundary conditions, in which case the Hamil- 
tonian is denoted by — 'Hr(c7) and reads the same as (2.1) except that the second term on 
the RHS has been dropped. 

The edges e = {x,y) come with a natural orientation: we will say that e is positively 
oriented if the vector [y — x) ■ > for all / G {1, . . . ,d}. Given a vector valued function 
/ : Z** — ?► R*^ we can associate a discrete vector field ( a function on edges) via 
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whenever e = {x,y) is positively oriented. In a region R and given f : R ^ M.^, we 
introduce the two discrete Laplacian operators by A^, (Dirchlet and Neumann) by 

- ■ fy: = fy where / is extended to be on 3°R, 

t/eRU3"R 

~ ^R ' fx = fx — fy 

x,y&R 

As a general notation, if / : R — > R*" is any (vector valued) function. 



For any A > set 



xeR 

fx=fx-f{R). 

r.^° = &'' = -(-AnA)-i..,, 

^ — ' — e(— Ar+A) • a-f. 



- &x,R 

Calculational Convention: In our estimates, we will use the constants C, c to denote 
generic universal constants whose values change from line to line. At various points in 
the work we use O(-) notation: An expression g G 0(/) for some other expression / if 
there exists a universal constant C > so that 

1^1 < c/. 

We will also say 

8<f 

if there is a universal constant C > so that g < Cf. 



3. Course-Graining, Contours, and Disorder Types 
We introduce two length scales £, L G {2'^)kGK where 

log2(L) = riog2(^-^|loge|4)]. 

The scales £ and L introduce coarse-grainings of R'^ and, by taking intersections, of 
Z'^. For any Lq G IN, we shall say that a block Qlo{^) is measurable relative to the 
scale Lo G N if r G LqZ'^. Denote this standard collection of blocks by Qlq ^rid set 
Af := ^o/igj— 32, ... 32, Y. For a fixed Lo-measurable box Qigir) and any ;^ G A/" let 
Qv := Qlo('- + ?7)- Set = {Q^ : Q G Qlo and G N}. Finally, Ql = {Q : 

Q- (^o/2,.../o/2) G QloI- 

For any set A C Z** we can associate a subset A C R*^ as the union of closed boxes 
of side length 1 centered at the elements of A. We shall say that A is connected if A is. 
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Note this is NOT the same as connectivity in terms of the graph structure on Z'^. We 
will refer to connectivity in the latter sense as graph connectivity. decomposes into 
one infinite cormected component Ext(A) and a number of finite connected components 
(Intj (A))j|Lj. The union of finite components is denoted by Int(A) := Uj!L^Int, (A). Let 
us denote the Lo-enlargement of a set A C K** by 

^Lo(^) = U{Q^^(,.)£(,.measurable:dist(QiQ(r),i)<^}Qi'o(^) (3-2) 

where dist(QLo()'), A) is the Hausdorff Distance between sets in R'^ in the £00 metric. 
We also introduce (^) / Irit, ( A), etc. by taking intersection of each defined set with Z'*. 
Most commonly, we will use this notation with Lq = L in which case we suppress the 
subscript, writing ^(A) The closed hull of a set A C is defined byc(A) := ^^(A)U 
Int(A) (with Lo = L). 

Given a spin configuration u E S we introduce the following phase variables: 



• Let 

for all r' such that dist(z, r') < 5i. 
otherwise. 

We expect (see § 8) that, typically, low energy configurations have 

the extra | log e| is taken here for convenience and is not optimal. In two dimen- 
sions, we must be more careful. 

• Let 

' 1 if (T{Qi{r')) • ei e [1 - ^, 1] whenever dist(z, r') < 5i, 
-1 if cr{Q({r')) ■ Ci E [-1, -1 + ^] whenever dist(z,r') < 5i, 
otherwise. 



= { 



• Let 



We will often suppress sub/ superscripts and the argument a from these functions, writ- 
ing i/jW,t/;(^),t/;. 

These phase variables are extended to the larger length scale L as follows. For coarse- 
graining purposes we consider only boxes Ql{^) which are L-measurable. If z G Qiij) 
with Ql{j) an L-measurable box, we set 



1 if t/;J'VJ°^ = 1 for all y G Qi^{r'), where Qi[r') G Qi 
and \\r — r' 



CO _^ 

= -1 for all y G Qdr'), where Qiir') G Ql 
and if II?' — / II CO < 2L, 



-1 iff?Vf 



otherwise. 



RANDOM FIELD INDUCED ORDER 
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By definition, Y is an L-measurable piece- wise constant function on A. 
Given C7 G 5, it is a standard argument that the connected components of 

R+ = {z G Z'* : Yf = 1} 

and 

= {z G Z'* : Yf = -1} 
are separated by connected subsets of 

R° = {z G Z"^ : Yf = 0}. 

Definition 3.1 An abstract contour T is defined to be the pair (sp(r), ipi) where sp(r) C Z^ 
is connected and L-measurable and ipriz) is a { — -valued function on sp(r) which gives 
the values of the phase specification on F. In the definition surrounding (3.2), whenever the 
set A in question happens to he sp(r) we will write 5{T),c{T), etc. Let Np = \5{I)\/L'^, 
= \3iT)\/i'^ so that Nf,N^ G N. 

Remark 1. (Other Choices of Boundary Conditions) Contours as we have defined here 
are not adequate for dealing with boundary conditions which are not ±ei because spins 
should be close to a chosen boundary condition in boundary boxes. One way to deal 
with this to redefine what it means to be good or bad according to how close a box is 
to the boundary of Ajv- In any case, the analysis of § 6 must be modified for boundary 
contours since we are not free to change the boundary conditions. 

Let us introduce the notation 3ext{T) = S{T) D Ext(r) and Si^{T) = S{T) n Int,(r). 
These sets are evidently disjoint. By definition of F, each of these sets is connected. 

Given an abstract contour Y we shall say that F is a concrete contour, or just a contour, 
for cr if sp(r) is a maximal connected subset of R^{o') and ^zicr) = ^t{z) on sp(F). It is 
possible that an abstract contour is never be realized concretely. 

We shall denote by 

X(r) = {cr G Q : F is a contour for cr} 

:= {a : sp(F) is a maximal connected subset of R^{cr) and ^z{<^) = ^r(z) on sp(r)}. 

(3.3) 

Put another way, we shall say that F is a contour for crifcrGX(F). 

From our definitions, specifying that F is a contour of cr lets us recover the values of 
Yz{cr),ipz{cr) on S{Y) just from the function tpriz) (see [25] for details). This convenient 
property allows us analyze systems of compatible contours without worrying about the 
microscopic spin configuration far away from the contour. Two contours Fi, F2 are said 
to be compatible if ^^(Fi) n sp(F2) = and i/^n = o^i the domain of intersection of 
5{T^),3{T2). 

So far we have considered contours at the level of spin configurations. We would 
like to be able to show that under certain finite volume Gibbs measures, a contour costs 
e^l\^\, the constant q being made large by appropriate choice of j6, ^ and e. In general, 
such an estimate will never be true uniformly in the presence of randomness. However, 



14 



N. CRAWFORD 



it turns out that large fluctuations of the fields, to be quantified in the next subsection, 
§ 3.1, and in § 11, don't effect the extraction of energy from most contours. 

Thus we turn to the interplay between spin configurations and randomness. We in- 
troduce (more) phase variables associated to the randomness a; G O. For this we need 
some more definitions. 



3.1 Disorder. Let Lq be a fixed scale of order e . Later we specialize to three scales: 
^/2,i and L as in (3.1). We will consider the behavior of the fields gQ for both Dirich- 
let and Neumann boundary conditions. The superscript distinguishing the boundary 
conditions is suppressed. 

First we need a coarse estimate on the behavior of the fields gQ. For a given cube Q 
we briefly introduce the potential 

e3x 

which plays an important role in §§ 5 and 6. We take A G [0, 1) fixed. Let A,B > Ohe 
fixed constants. Consider the event 

AQ{r)= \co : r-^ ^^Q,y ^ all x G Q s.t. distco(x, 9°Q) > ^ > • 



\\y-x\U<r J 
A box Q G Q^^, will be called nice if for all A < Lq ^: 

The event Aq^^ (r) occurs with r = log^° Lq. (3.4) 

||g^J|eo<eL-^|loge|30 (3.5) 

l|VgJj|co<e|loge|30 (3.6) 

^e'< l|V^Jjli/|Q,| <Be2 (3.7) 

lklU,Q„ < llogepo. (3.8) 

kQj^Q^< |loge|3°. (3.9) 



A box Qh^ir) G Qlq/ will be called good if Q,^ is nice for all rj E Af. A box Q G Q'^ is it 
called good if it intersects a good box in Qig. We introduce the function on Q'^ by 




1 if Q is good, 
otherwise. 



Appropriate choices of A, B > in (3.4) and (3.7) will be made below. We make 
appropriate choices below. The factor | log ep'^ is so that boxes which violate the bounds 
are exceedingly rare if e is small. We have not attempted to optimize this part of the 
proof. All of these requirements must be modified in two dimensions. Also, the bound 
(3.5) may be, in principle, tightened as the dimension increases. 
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Let y be an Lo-measurable connected and bounded set. Suppose that F, G are real 
valued functions on the collection of cubes in Y which are Lo-measurable and define 



[f;G]y= X] f(Q)G(Q). 

We will say that Y is good if 



QcY 



[1-Evl]y= X] l-ELo(Q)<|loge|-55^^° (3.10) 



QcY 



where N.^" - 



For notational convenience we introduce the functions on Lq measurable blocks 

12 . l|Vgi„2 



fAlQio) = max ' , (Qlo) = max ■ 



K(Qlo) = max||a||oo,Q,,. 



An Lo-measurable connected and bounded set Y will be called regular if the following 
estimates hold for A G {0, Lq^^ log Lg}: 



pV. 1 

A '-^{f7>e2|loge|} 



^ < emy° (3.11) 



^'^'"^{FA>e^[A-2ALo]|loge|} 



<|loge|2Ny'' (3.12) 

Y 

[K';lR>|iog<rpo]y < lloger'^N^" (3.13) 

E E \<Qri)\ < ^o' los'° ^oN^° (3.14) 

Lq— measurable 

Definition 3.2 (Clean Regions) An L-measurable connected and hounded set Y will he called 
clean if it is good and ifS{Y) is regular at the scales Lq G {^/2, i, L}. Otherwise, call the region 
dirty. 

Remark 2. Not all conditions are used at all scales. Its just faster to formulate our re- 
quirements in a uniform way. 

Remark 3. (Two Dimensions) These bounds are not sufficient if d = 2. In this case, we 
need to control ||g'^||oo, || Vg'^||oo more carefully and also add requirements regarding the 
size of 

E fA^(Q)' E fA(Q). 

QCY QCY 

The following proposition is fairly straightforward. The parts of it we prove are post- 
poned until § 11. 
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Proposition 3.3 Let Y be connected L-measurable and bounded. There exist A,B > so that 
the following holds. We can find cq G (0, 1) and a constant c > so that ifeE (0, eo) 

P (Y is dirty) < g-cliogepN^. (3.I5) 

Let 

^ = {y C Z'* : y is connected L-measurable and is not clean}, 

D = UYe^c(y), andDA = DnA, 

where A is a connected, bounded L-measurable subset of "ZJ^ and 

For any x ^ "Z? ,\e\. Q{x) denote the L-measurable box containing x and define 

A{x) = {a; : Q(x) is in c{Y) for some Y which is not clean} . 

We relegate the proof to § 11. 

Lemma 3.4 Choose A, B so that the estimates in Proposition 3.3 hold. There exists eo > and 
c > so that ife E (0,eo) 

Var{r)A) <e-'^°s^^\A\ (3.16) 
Further, for any sequence t ^'^/ ^ ^-S- there exists Nq{lo) > so that 



A 



N 



for all N > No{co). 



4. Statement of the Peierls Estimate and Proof of Theorem 1.1 
Given a contour T, T* will denote the (O, P) event 

r* = {o^ : ^(r) is clean and c{Y) is not strictly contained in D} (4.1) 
If r* occurs, r will be called a * -clean contour for short. Given a spin configuration 

x(r^, . . . , r*,, Tm+\r ■ ■ ■ / ^m+n) = n,x(r/) 

where (Fi, . . . , F,„) satisfy the event defined in (4.1) and (F„,+i, . . . , F„j+„) do not. Oth- 
erwise we define the right hand side to be the empty set. As a variation of standard 
definitions, let us say that 

(F^, . . . , T*^, Tm+l, ■ ■ ■ , Tm+n, 

are *-compatible if 

Fm+l/ • • • / Y m+n) 7^ ^ ■ 

Lemma 4.1 (The Peierls Estimate) There exist 5 E (0, 1) so that ifO<^<S,we can find 
£o = £o(^) so that the following holds. Ife < eo there exists fie so that for fi > fi^ we have the 
estimate: 
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Let N be fixed and consider the event JL{T\, . . . , F*,, r„)+i, . . . , r,„+„). Then 

f/X^(x(n,...,r:„r„,+i r,„+„)) < e-^^^^^\^'\ 

where 

g = CfV|loge|-24 (4.2) 

for some constant C > 0. 

This lemma is proved in several steps below: §§ 5, 6, 7 and 8. Theorem 1.1 is then 
completed via the following Peierls contour counting argument. 

Proof of Theorem 1.1. The proof of Theorem 1.1 is an application of Lemmas 3.4 and 4.1. 

Let A, B > be as in Lemma 3.4. Let x is taken so that Qiix) CiD = 0. Note that this 
does not depend on A^- If t in the Van Hove sense then eventually Qi {x) C A^. 
Consider the event 

{c7:Yq^(,)(c7)^1}c5a,. 
By definition of xp, Y and the ei boundary condition, there exists a largest contour T so 
that Ql{x) C c(r). Moreover, since Ql(^) H D = 0, F must be clean. Decomposing 
V^Qi^(x) 7^ 1} irito disjoint subsets according to this largest contour we have: 

T,w *-compatibIe: Qi{x)Cc{T) 

By Lemma 4.1, 

as long as q > 21og(2flo)- The theorem now follows easily. □ 

5. Approximate Ground States in the Bulk of a Clean Contour 
Let r be fixed and assume cr G X(r). Let 

s{Y) = Sl,^{y) n Am, 7^f,^(^) = {Qe Q<,^ Q n s{y) + 0}. 

Given a region K C A^ let ext denote the boundary condition which is set to e\ on 
3°R n A^ and is free otherwise. With this definition set 

Eo(Q) = max-'HQ(cr|ext) 
The goal of the present section is the following Lemma: 

Lemma 5.1 (Bulk Ground States) Suppose that T is a clean contour. There exists eo > 
and a constant C > so that if < e < cq the following holds. On S{T), there is a spin 
configuration cr^^^^ which satisfies 

< E UQ)+^siT){<^"-'^) < e'lloger^^lrl 
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Further, 

I - 4'^ ■ ei < Velos'' e 

for all y G 5{T) and 

ll^"('^-^2|U,3,,-(r)<e|logep. 

In the remainder of this section we give the basic construction. In § 10.1 we estimate 
the errors made in the construction. The construction proceeds by stitching together 
approximate ground states on the microscopic blocks Q G TZt/JJ^) There are two cases 
to consider. 

Case 1: Sf/^(Q) = 0. We then set ay = ei for all y G Q. This is unlikely to be (even 
nearly) optimal in Q but because we are working in a regular contour this error will not 
have a large contribution when calculating the energy difference between the configu- 
ration we are constructing now and any a G X(r) (this follows from a combination of 
Proposition 8.1 and (3.11) and (3.14)). 



Case 2: Ef/^(Q) = 1. Provisionally we define, for all y G Q, the angle 
and the spin configuration 

CTQ^y = (cos(0y),sin(0y)) . (5.1) 
In order to stitch configurations in neighboring boxes we interpolate the function dy in 
Q so that the angle is small near the boundary of a box. Let 

dist(x,3°Q) , 

if X G Q and 

at}^^ = (cos(Ty 0j/),sin(Ty 0y)) if y G Q. 

The reason for this slightly complicated definition is that it guarantees \Ty Qy\ < 2£^^ for 
all yGB^Q. 

6. Surgery at the Boundary of a Contour 

Suppose a contour F is given and that a G X(r). These will be fixed throughout this 
section. Let 

c(r) = ^(r) n sp(r)^ n Am, c±(r) = {z g e:(r) : Y,{a) = ±1} 



i.e. £(r) is a thickened version of the boundary of F. Consider also the 'middle' portion 
of (2:(F): 

m{T) = {z G £(F) : dist(z, [9°e:(F)\A^^]) > | - 100}, M{T) = {Q G Qf,^ : QnM(F) ^ 0}, 
M(F) = Uq^m{t)Q. M±(F) =M(F)n£±(F). 
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We will often drop reference to F in these sets below as F is fixed throughout. 
The main result of this section is the following. 

Lemma 6.1 There exists eo > so that the following holds for all e G (0, eo)- J^^i T bea clean 
contour. 

Let a G X(F). From cr, we can construct a spin configuration a'^ with the following properties. 

(1) (7-,. = z/dist(x,e:(F)) > f . 

(2) Ifx G d"M^{T) then ||c7-^^ei||2 < e|loge|. 

(3) For any choice of boundary condition Ci compatible with F and for which sp(F) C A^, 

- -HA^^'^ki) > -HaM^i) - ce'l iog^r''|r|. (6.1) 

The proof of this lemma consists of four separate modifications of a. We will first 
present the essential steps and Lemmas used in the construction. Proofs will be given in 
later sections. 

6.1 Modification 1. Let 21^ denote the smallest set containing £^ so that 

s^n{a^ ■ ei) = ±1 for all :r G 3°2l^. 

The following is a consequence of Lemma 9.2. Recall that ^ defines the cutoff for whether 
block averages of spins are good or bad. 

Proposition 6.2 There exist > and eo(^o) so that ife G (0,eo) cind ^ G (0,^o)/ 

2i± u a^a^ c {x : dist(x, < l} 

On the basis of this fact it is convenient to make a first modification of a, into crj, as 
follows: If X G 21^ and sgn(£7;c ■ £i) = Tl/ set cr^ to be the reflection of across the 62 
axis. Otherwise set al = Cx. Let us record, without proof, some simple properties of this 
transformed configuration. 

Lemma 6.3 Given cr G X(F), let be given as above. 

(1) G X(F). 

(2) Ifx G e:±, sgn(t7-i . ^ ^1 if^^ • ei 0. 
(3) 

k'(Qf)-ei| > \a{Qe)-ei\ 

whenever Qf c C 

(4) IfY is compatible with the boundary condition ei, 

--HA^cr'lei) > -UaM^i)- 
For any region R C An 

£R{cr') < SR{cr). 

(5) The number of spin configurations ct" G X(F) which map to the same is bounded by 

Irl 

for some universal constant c^ > 0. 
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6.2 Modification 2. From now on, we will work with a^. Ideally, we want to use 
Lemma 8.2 to show that by replacing cr^ by the optimizer of — 'He:(c^|c^) in £ the re- 
sulting configuration behaves as claimed in Lemma 6.1. However, because the random 
field may behave rather poorly in ( relatively small) subsets of C this cannot be done 
directly. 
Set 



^ - UQeQj^:Qce:andSt(Q)=oQ 

& = {x: dist(x, D) < 5L}, = . 

S± = {x G : dist(x,9°i^±) > L'}. 
In the following, we work with the regions 



'nc^(r). 



(6.2) 
(6.3) 
(6.4) 



To prepare the way for Lemma 8.2, we first we will brutally change in . Define 

^ Al forx G &^\^t , 



16dist(i-,2)± ) 



ioxxeT>t , 



■he' 



otherwise. 

Representing (7^ = (cos(0v)/Sin(0;(.)) with 0^. G 
(resp. if X G ^ {I)) let 



[-f,f] (resp. G [f,f ])if X G e:+ 



if X ^ 



(cos(T;,0;,),sin(T^0O) ifa:G^+, 
^{cos{Tx[Qx - tt] + 7r),sin(Tx[0T - tt] + tt)) if x G 

Lemma 6.4 T/zere is so that ife G (0, eo) and ifT is a clean contour, 



< 



loge 
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6.3 Modification 3. Set 

Our next goal is to modify cr^ in 'lo^ in such a way that the resulting configuration has 
small projection onto the vertical (ea) axis for all x G 9°M without too much cost in 
energy. Note this aim was achieved already within 2)^ . 

To be concrete with our calculations, we restrict attention to ; the region 

is treated similarly. Let g'^'^.+ = [— A + A]^^a with dirichlet boundary con- 
ditions on 3°^+ and with 

A = L-2log^L. (6.5) 
Because of Modification 1, for x G U 9°"^+, we can express in angular variables 
as (7^ = (cos(0x)/Sin(0j)) with 9x G [— f , f ]• Consider change of variables 

(px = ex- cos(0;,)4'J+ for X G 
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We show later that ||oo ^ using Lemma 9.4 and the definition of '^+). Thus 
the transformation is nonsingular and preserves the half-space (7 ■ ei > in spin space. 
Let 

^(r) = {x G e:+ : dist(x,9°e:+ n An) > Vs} 

and let f be the smallest subset of 1? containing ^ and so that 
for all X G 9°f. 

For e small enough, our definition of £, Lemma 9.2 and the fact that ||oo ^ ^In 
together imply the following. This is the second point at which we use low dimension- 
ality of the lattice 1? . 

Proposition 6.5 There exist > and eo(^o) so that ife G (0,eo) ^ G (0,^o)/ 

f C {x G : dist(x,3''(!:n An) > V9} . 

Set g+ := f H After the change of variables in according to Lemma 8.2 it is 
sufficient to optimize (after restricting /C^+ to 0+) 

^s+if\4>)= Li cos{Ve(p) + J m^cos^icpx) (6.6) 

where = Ei/~x [S'f^'? y ~ i^i^o^- (6-6), we neglected the boundary terms in (8.3) and 

also the Error^j^l term. The former plays no role in the optimization and the latter will 
be bounded carefully in the proof of Lemma 6.6. The advantage of (6.6) over the original 
Hamiltonian is that stationary points of /Cg+ (• j^) behave as solutions to discrete elliptic 
PDE's with (random) mass. 

Let !? be the maximizer of /Cg+ (^|^) on g+ and let 



if X G 0+ 
(/)xifx G £+\0+. 



Inverting the change of variables, define ^+ by 

<^;. = 0+-cos(0.+ )^Jf^^ 

We construct the configuration <^ close to tt in a similar way on 0^ (fluctuating around 
the angle tt). Modification 3 of the original spin configuration is defined by 



(cos(0+),sin(0+)) forx G '^+, 
(cos(4>^),sin(0^)) forx G 
ai otherwise. 



Lemma 6.6 There exists £0 so that ife<£o 

-JiAA^^h) > -nAA^\ei)-Ce^\\og£\-^^\T\. (6.8) 
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Moreover, 



= (p- for X G An\0^ U g , 
1^.1 -4'?J ^ exp(-clog4e) ffdist(x,M+) < V5, 
flnd similarly for <t>^ . 

6.4 Modification 4. The configuration is almost what we want. We make a final 
modification into (7*^ so that the resulting configuration is close to ±ei on the appropriate 
components of 9°M. 
By construction 

OX := I X G (J:(r) : dist(x, d°€, n A^) CSt Ug+UST Ug" 

[ 5 J 12 12 

and, obviously, M is contained in Let (c,);g[j] denote the connected components of 01 
and letM, = c,nM(r). 

To achieve the stated goal we interpolate in a way similar to that employed in § 5. For 
X G ^± let 

1 otherwise. 

Let the angle dx be defined by = (cos(0;c),sin(0x))- If c, C £^(r) we may take 
dx G [-f/f] and then set 

ai^x = {cos{Tx9x),sm{Tx9x)). 
If C; C (t"(r) we may take for Ox G [f , ^] and then define 

= (cos(Tx[0x - tt] + 7r),sin(T;,[0x - tt] + n). 

Finally we set 

_ /c^/>if ^ e c,n(5, 

1 otherwise. 



Lemma 6.7 There exists eo so f/jfli if e < Cq and if! is a clean contour. 

nA^{(r'\ei)-n^,{a^\e^)\ < e^] ioge|-40|r|. (6.9) 
Moreover, ifxE d°M^{T) then ||cr^ T S1II2 < ^| loge|. 
Let us finish this section with the observation: 

Proof of Lemma 6.1. The statement of Lemma 6.1 follows from the construction, collecting 
the estimates from (6.8) and (6.9). □ 



7. Gluing Configurations and the Proof of the Peierls Estimate 
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Given a G X(r), the next step in our construction is to attach C7^(^' as constructed in 
§ 5 to a modification of the configuration a'^'^^'^ as constructed from a in Lemma 6.1. 

We shall say that T is a ± contour if Y2(r) = ±1 on SextiJ) (recall that §extij) is the 
exterior boundary component of F as defined below definition 3.1). To obtain the desired 
configuration onc(r) = ^(r) U Int(r), we distinguish whether F is a ± contour. 

Restricting the spin configuration ct"'^'^) to A]v\(5(F), we produce a new spin configu- 
ration a* on A^\5{T) as follows. If F is a + contour, then on each interior component 
with Y^'^'lc/ ,-r.-. = —1 we reflect all spins across the e^-axis. For each of the remaining 

components of A-f^\5{T) we keep the spin configuration fixed. 
We define 

If instead F is a — contour, we reflect spins on interior components with Y*^'') = 1 and use 
the reflection of ay across the 62 axis in place of ay , calling the resulting configuration 
Sp. The proof of the next lemma appears in § 10.3. 

Lemma 7.1 There exists £0 G (0,1) so that for all e G (0, eo) the following holds. Suppose 
that F is a clean contour and that a G X(F). Then 

-■HA,(S±((7)|ei)+7^A,(c7|ei) > ^V|loge|-24|F| (7.1) 



7.1 The Peierls Estimate. We are now ready to derive a Peierls estimate for this sys- 
tem. By definition, if (7 G X(F), the restriction of ^z{o') to each of the components 
Sext{T),3l{T) is constant. Leti^±(F) = -R^\^^(F). 

In order to deal with contours adjacent to let 5n{T) = S{T) n A^. Given that F is 
a clean + contour, we say that a spin configuration cr<)-^(r)' G ^^^{yy is compatible with F 

if '^^N(r)^lA5v = ^1 ^^^^^ 



K{T) [^M,(T)) = Mz) for z G sp(F) and Y,((7;^(r)) = ±1 forz G K±(F)j ^ 

Note here that implicitly takes as an argument the extended configuration (c^'^^i^y)' ^^n{ty) 
although we will suppress this detail below. 
For any such C7^^(r)c/ let 



^^rT ("P-KiT)) = '^r(z) forz G sp(F) and Y,(c7-;^(r)) = ±1 f^^ z G R±(F) 



l^s^ij)'"''' (Y.((7;^(r)) = lforzGMr) 

Here o'lj^^ijy denotes the boundary condition obtained from Lemma 6.3 and ^+ ■ c^^^^p^ 

is obtained by making a global reflection of a^ on each component of R n [c(F)\^(F)]. 
Let 

||W(F;-)||= sup W(F;^,^(r)0. 

Ws-^{TY compatible} 
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Notions for — contours are defined similarly with the provisos for + and — reversed. 

Lemma 7.2 There exist 6, cq, /3o, C > so that if e < Cq, fi > fig > fio, < ^ < S and T is a 
clean contour 

||W(r;-)|| < 

where 

^? = c/3^viioger24. 

Proof. Assume for concreteness that F is a + contour. Given a configuration o'Si^{ty com- 
patible with r let c7-;^(r) ^ 5^^(r) so that tj' := {(r'^^ijy(r5^{Ty) e X{T). Let S+K) 
denote the output of our construction from just above Lemma 7.1. Dropping the su- 
per/ subscripts, 

where is the first modification of a' as defined in § 6. Note that S((7')|^j^(r) is an 
element of the event in the denominator of W{T, o-^^i^ry)- 
Let 

F^, = {ae Ss^^T) ■■ \\(rx - S{(/),\\2 < Vx G ^M(r)}. 

Then 

I^Zirf'' {"^^KiT)) = 1 for z G S^{T)) > ^fj^^^^' (F) . 
Using (3.13), we have 

xG<5(r) 

Using this bound, it follows from the definition of F^i that 

By Lemma 7.1, we therefore have 

W{T,crs^^y)<e^-f'f^'^+s^'^^\^\ 

where 

/(e) = ci?viioger^4 

accounts for internal energy difference between the numerator and denominator and 

g{e) = C2|loge| 

accounts for the entropy difference. The lemma follows immediately. □ 

Proof of Lemma 4.1. Let N be fixed and consider the event X{Tl, . . . , F|^, Fm+i, • • • , ^m+n) 
where the contours are assumed compatible with the boundary condition ei . Then we 
claim 

m 

Vm+i F,„+„)) < nCW(r/;-)ll 

where is the constant from Lemma 6.3. Once this is justified, the Lemma is proved by 
application of Lemma 7.2. 
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The proof of this claim proceeds by induction on m. Interpreting an empty product 
as 1, the case m = there is nothing to prove, so we proceed to the induction step. 
Suppose the claim is true for any *-compatible system with m = k clean contours and 
n dirty contours. Given a set {FJ, . . . , F,*^^, Fj;+2/ • • • / ^k+n+i)} of *-compatible contours 
and reordering as necessary, we may assume 

c(F^+i)nutisp{F,} = 0. 

Assume for concreteness that Fj.^i is a + contour. The argument in the case of a — 
contour proceeds in a similar manner. Let {f i, . . . , f r} denote the set of contours among 
{^m+i,- ■ -f^m+n} with ^^(f,) C Int(Fjt+i) and . . .,F'^_^} denote the rest. Then we 
have 

x(ri, . . .,F^^;^,F^+2/- • •/Fj:+„+i) = x(F5|;,. . .,f|,f'i,. . . ,y'„_,.) nx(F,*_|_;^,ri, . . . ,Tr) 

Using the DLR equations we have 



Note here that 

^^^in+v^i r,)/cM(r,+i) 

may be depending on t^c^(r,+i)- 

If r = (sp(r), ipr), let -r = (sp(r), -ipr) and let 



—ff otherwise. 



Recall the map a i-^ from Lemma 6.3. It is straightforward to show that 

/I 

\-^^(r,:+iTi r,)/cM(r,+i) 



[pi 

Because the map is at most -to-1 for some universal constant > 0, (see Lemma 6.3), 

M^t n,T[ r;,_,);ix(Tr,^j(ri) t,^^,{t,)))c^(y,^,) 



N 



<ci[^+^i||w(r,+i;.)||?^^,(x(r5[,...,r|,r; rU,Tr,,,(fi) Tr.Jfr))) (7.4) 

and the induction step is proved. □ 



8. Energetic Estimates 
This section is the engine room for the entire paper, see in particular Lemma 8.2. 
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Proposition 8.1 Let us suppose that L' G N is fixed. Then for any A > and any cube Q of 
sidelength L' 



^E^e>-cr, + 0(e|4Q)||Q|). (8.1) 



Similarly, for any finite region R d'ZI^ and any A > 



1 ^ r„ An^^2 



E V.gR ) + A E ^R:. ^2 • + E (8.2) 



enR/£ 



jGR 



x~j/,.tg3'R 



Proof. Consider the free boundary condition case A = 0. Solving the equation —A^g = a. 
in Q, we can write 

-^q(^) = -Uq{^) + E Veg^e2-Ve(T + 0{e\a{Q)\\Q\). 

ecQ 

Completing the square gives the expression first expression. The remaining cases are 
similar. □ 



For notational convenience, let us write cr^ = {cos{6x),sm{6x)) with 9x only defined 
modulo 2n. For any pair of nearest neighbors x,y, let [Veff)]!^,^ denote the minimum 
of [Vg^)]'^ over all pairs of angles in the equivalence classes which give the vectors 
(cos(0x),sin(0,.)). (cos(0j/),sin(0y)). We have 

((7,-^7,)2x[V,0)]L. 

Given a region R C Z'^ which is finite and any x ^ R, suppose that we represent a 
spin Cx = (cos(0j), sin(0x)). We will find it revealing to make the change of variables 

O'x =0.-g'x 

where g'x = cos{ex)g^l^ 

org; = cos(0,)g^;^ 

depending on the context. Note that despite the fact that 6x is only defined modulo 2tz 
this transformation is unambiguous. Moreover, it is nonsingular as long as \gx\ < 1. 
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d 



Lemma 8.2 Let us suppose that R C is fixed and hounded. Then for any A > 0, 
-^rM'7)= E [cos{V,e')-l] + -Vm,cos^{e',) 

er\R^O xeR 

+ E g^'''sm(0,) + i E cos'{ey)[g'y-''-g'.'''? + Erro^/'{a) (8.3) 

i/e3"R y&dOR 
x~y,xed'R x^y,xGd'R 

where 

mx=Us';'''-s'x'''? 

yr^x 

and 

Error'^'' < M\g'^''h\\e2 ■ ^Ur + (|| Vg^'°|U + l|gnU)(ll V^||1«u3"R + II 111)- 

(8.4) 

Similarly, for free boundary conditions we have 

- ^Ria) = E [cos(V.0')) - 1] + 7 E m.cos'ie',) + Error''/' (8.5) 

eCR xGR 



{x,y)cR 



and 



Error'/" < 

eM\\a- e2h,K + A|!^^'^||2|k2 ■ ct\\2,r + (|| V^^'~|U + Wg/'UMVcrWlK + WVg/'Wl). 

(8.6) 

Proof of Lemma 8.2. In the proof we use the notation fe = ^'/' where e = {x,y) for any 
(vector valued) function / on vertices. 

We prove only the first estimate. Let g'^ = cos{6x)gx'^ and note that g'^ is independent 
of the branch chosen in the definition of Ox- 

First, for any edge e with e n R 7^ choose 6x and Oy which achieve [Ve0)]27r- This 
choice is, in general, edge dependent. Letting 0'^ = 6x — g'x^ 

cos(V,0) =cos(V,0')-sin(V,0')V.^'-^cos(V.0')(V,^')' + O([V.gf ). (8.7) 
For a positively oriented edge e = (y, x), we will use the formulas 



9x = 0e + -VeO 



28 N. CRAWFORD 

Using them, we may write 



V,g' = COs(0,)V,g^'^ + O(|[V,J]2,r||^^'°L+ [V,^^'^]2) (8.8) 
Combining (8.7) and (8.8) and summing over edges 

X: cos(V,0)= cos(V,0')-cos(0,)V,0V,g^'° + ^^4^(V,/'°)2 

enR/0 enR^0 

+ 0((||V^^'^||oo+||g^'°||oo)(||Vc.||2j,uaoR + l|Vg^'°||i)). (8.9) 
Now consider the contribution to — "Hr coming from the random field. We have 
e ^ a, sin(0,) = ^ (-A^ + A) ■ g^,'^ sm{9,). 

xeR xgR 

By Cauchy-Schwarz, we have 

I J2 A^^'^sin(0,)| < A||g^'^||2,R|| sin(0)||2,R. (8.10) 

xeR 

Making a summation-by-parts, 

X:-AD-g.^°sin(0,)= E V,/'^V,sin(0)+ g.^°sin(0,). (8.11) 

xGR er\R^0 yeSOR 

3r~j/,xe3'R 

Focusing on the first term in the RHS of (8.11), for each edge e = {x,y) we choose the 
same branches of 9x, Qy as in (8.9) and expand around Og: 

E V,/'°V,sin(0)= E cos(0,)V,g^'°V,0 + O(|V,/'^|[V,0]L) (8.12) 

er\R^2i er]R^0 

If - edge by edge - we combine the RHS of (8.9), (8.11) and (8.12) we obtain (8.3) with 
second term replaced by 



(8.3) is obtained by first replacing — with ^ 
bounded by the error already amassed and noting 

,2 



cos2(9,) , cos2(ey) 



at the cost of a term 



E ^^^m. = i E cos^(^^;)'«. + O E I^IJV.g]^ . 

xeR xgR VeGR / 



□ 



9. Preparatory Lemmas for § 10.2 



9.1 Absence of Defects: Proof of Propositions 6.2 and 6.5. In this section we consider 
the extent to which the Dirichlet energy Eq^^ {a) can be used to control 'smoothness' of 
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spin configurations in cubes Qig (these need not be Lo-measurable. We consider this 
question in Z'^ for d > 1. 

Definition 9.1 Let }i G (0, 1) and let < S < }i be fixed. We shall say that a G Sq^^ has a 
defect (iv.r.t. 5, }i) ifcriQig) ■ ei > 1 — 3 but there exists a graph connected subset R C Qlq so 
that diam(R) > ^ and dist(R, QLoi^Y) > t ^""^ that for all x E R, ■ ei <\ — }i and 

We wish to argue that defects cannot occur for low (Dirichlet) energy configurations 
for d G {2,3}. The following lemma depends on low dimensionality. In § 6 it allows 
us, at the boundary of clean contours, to restrict the Hamitonian to regions of the phase 
space on which it is (morally) convex. 

Recall I ~ e^^l logej^'*. A similar statement can be pushed through in the two di- 
mensional case. 

Lemma 9.2 (Energetic Cost of Defects) Let d = 3 and fix < S < ji < 1. There exists 
£ G (0, 1) so that if£Q{a) < logel^'', cr does not have defects in Q^. 

This lemma is an easy consequence of the following. Let B/ denote a fixed £^ ball of 
radius / in Z**. 

Lemma 9.3 (Point Defects in ^-Dimensions) Let d > 1. If ji €1 (0,1) and S G {0,}i) then 
for all I sufficiently large, whenever 0'{Bi) ■ ei > 1 — 3 and ctq ■ ei < \ — }i 

(}ifd = I, 
W>\^,ifd = l, 
U ifd > 3, 

where the implicit constant depends on 5, }i and dimension. 

Proof of Lemma 9.2 given Lemma 9.3. Let us assume for convenience that Q( is centered at 

0, that is Q^ = {-|,...,f-l}'^. 

In let {^1,^2/ ^3} denote the usual basis of orthonormal, positively oriented imit 
lattice vectors. Let Hj{k) denote the hyperplane through kci perpendicular to e, and let 

Bi{k) = Hi{k)nQ(fork G {-^/2, • • • //2- 1}- We may decompose £7(Q^) = E^^^ (t{B 

1. e. as a sum over d — 1 dimensional hypercubes perpendicular to each e,. Setting 

fi{k) =(T{Bi{k)) we have 

r'j:\m)-a{Q,)f<r' lkx-(7(Q,)||iand 
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Inequality I follows from the Poincare inequality and II from our hypothesis on the 
Dirichlet energy By Markov's inequality and the defintion of I, for each / 

^-^\{k■.\\m-a{Qf)\\2>l^-^/2}\<y,,-s\\o^e\-\ (9.1) 
^-'\{^--£Bm^^) > llog^mi < |loge|-i. (9.2) 

Thus if Qi has a defect, we can find at least ^/s disjoint 2 dimensional subcubes Bi {k) so 
that 

fi{k)-ei>l-[S + F-i/2], (9.3) 
£B,{k)icr) < |loger^ (9.4) 
cT^c ■ ei < 1 - for some x G Bi{k) so that dist(x, d°Q() > ^U. (9.5) 

For each such x let hx be the P- ball of radius around x in i3, (A:). Then 

For e sufficiently small, this is in contradiction with Lemma 9.3. □ 

Proof of Lemma 9.3. Let us assume for convenience that B/ is centered at 0, i.e. B/ = {x G 
• W^Wi < The proof of the Lemma is easiest when d = 1: The condition c(B/) ■ Ci > 
1 — S implies that we can find a pair of vertices x < < y such that Cx ■ Ci, ay ■ ei > 
1 — [S + I'^^/i] and x,y E [—/,/]■ Optimizing the Dirichlet energy £ in [x,y] subject to 
these boundary conditions and the condition ao ■ ei < 1 — }i then gives 

For dimension d > 2 it is most convenient for us to replace dx by a continuous interpo- 
lation. For y G R let [y], {y} denote the integer and fractional parts of y respectively and 
extend these notations componentwise to vectors in K^. Let Bj^ = {x E M.'^ : [x] G Bj;}. 
We define a on B/_2 by 



ve{o,l}'i 



YI(V,{X,} + {1 - - 



/=1 



i.e. d'{x) is a multi-linear interpolation of c^"'. 

By definition, we can find a small ball B' (depending only on d, 5 and }i) around so 
that a{x) ■ei<l-[5 + ^i^-^)/^ if x G B'. In addition, 

\Bi-5\'^ f i7{x)-e,d'^x=\Bi\-^ J^(T-e,+0{yi) (9.6) 

•'^i-s xeB, 

with the second term on the RHS coming from boundary integrations. Taking I large 
enough, we may assume |B/|^^ Jg^ d'^xd'{x) ■ ei > 1 — [S + (h-^)/4]. Further, there exist 
universal constants Ci,C2 > so that 

ci E [V.^^]'< /- d'x\Vax\^<C2 E ["^ecr]^ 
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where Vo'^ is the gradient of in R"^. It is therefore enough to obtain bounds in the 
continuum setting with |cr;i;| = 1 relaxed to \d'x\ < 1. 

We assume d'x ■ Ci > Q and that cTx is spherically symmetric (although perhaps not of 
unit length) since the operation of taking spherical averages preserves (9.6) and can only 
lower the kinetic energy. 

Now we proceed in a manner similar to the one dimensional case. We can find a 
sphere 

S, = {x G R'^ : ||x||2 = r} 
with r < / so that 5-(||x||2) ■ ei > 1 - + i^-^/i]. 

Combined with the condition P'(x) ■ ei <\ — [5 + 3(/'-<5)/4] if x G B' and using capaci- 
tance estimates between balls in R*^, this implies 

for d = 2 



h, ' ' ~ 1 1 for d 



J for d > 3 

The implicit constants depend on the capacitance between a unit ball and one of radius 
two and also orv5,}i. □ 

Proof of Propositions 6.2 and 6.5. These are easy consequence of Lemma 9.2 □ 

9.2 Localization Estimates for Given a region R, in this section we show how to 

control the field g'^^ in terms of the family of fields {gQ^ ) q where the boxes Q are of the 
form Q = Q{r + rj) for some Q(r) G Qi and some rj E N := j^{—32 . . . ,32}-^. Here, as 
in§ 6, A = L-^log^L. 

Recall that in § 3, we defined the function R on Qi by R(Q) = max,^ II*I|c»,q,,- Recall 
that Q'l^ = {Q}] '■ Q ^ Ql arid rj G JV}. Here is a convenient, though not particularly 
optimal, bound which we use later: 

Lemma 9.4 Let Q G Q'^ so that QCiR 0. Suppose there exists S > so that x G Q n R 
and dist(x, 9''[QA^]) > A^zllogel"^. Then there exists c > (independent of S) so that 



QnR^0 



-c\log€\^ ^ ^-cAJ dist(x,Q 



In particular, i/supg^^^^^ K(Q) ^ ^ 5 > 1, 

Let us emphasize here that 3° [Ri AR2] denotes the outer boundary of Ri U R2\Ri H R2- 
In particular, it is possible for x G Ri H R2 to satisfy the hypotheses of this Lemma while 
being close to the boundary, as long as the part of the boundary it is close to is shared by 
bothRi,R2. 

Remark 4. One annoying problem if d = 2 is that this bound cannot be used as the 
restrictions on £, L are too tight to take ^ > 1. 
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Proof. The proof exploits a random walk representation and coupling argument. Set 

Ri = R,R2 = Q. Let Xj'' for i G {1,2} denote continuous time simple random walks 
started from x. In order for things to work out correctly, we take our walks to have 
exponential holding times of rate 2d. For each / G {1,2}, let T, denote the exit time of 
X(') from R, and let T\ be an independent exponential random variable which has rate 
A. We couple X^^, X^^ until the first time the coupled walk hits Ri AR2. Call this hitting 



time T3 (it can be infinite). Then 



iRi,x 







o.A,D _ A,D 
SRi,x Sr2,x 



axil{Ti < Ta,T3} -E.V ax2l{T2 < Ta,T3} 



I II 
Clearly 

II < R(Q') E [t21{t2 <Ta,T3}] 

for any Q' G Ql which intersects Q. We claim 

E[T21{T2<Ta,T3}] <A-V^|l°g^r 



(9.8) 



Note that on the event {t2 < T3}, \Xj^ — ^| > A^zllogel"^. Thus,ifr(A) = V2A^2|loge|^ 
and T,.(A) denotes the first exit time from the ball of radius r(A) around x, we can 
bound the expected value by 

E[t21{t2 < TA,T3}] < E[ta1{t,(a) < Ta}]. 

The inequality (9.8) now follows. 

Term I is estimated similarly. Let Tq be the hitting time of Q G Qi for X^^^ Then 

I < E R(Q) E,. [ti1{tq < Ti < Ta, T3}] . 

Reasoning as above, 

E, [ti1{tq < n < Ta,T3}] < A-ie-^|i°g-i'^ Ae-^^^'^i^'^'Q) 

from which the claim follows. Note that we are being quite generous here as there 
should be cancellations which reduce I and II. □ 



9.3 Estimates for Maximizers of /Cr. In this subsection we consider the behavior of 
optimizers of the functional 

^A(P\'P)= E cos(V,<p)-l + J X] "^-^-cos'l-p.) (9.9) 



RANDOM FIELD INDUCED ORDER 33 

where = J^y^xis^? y ~ Sv^]^- summarize, the properties we will use below are: 

^I3°Aa,U[S)+ n3»g+] = 0. 

'12 

Lemma 9.5 Suppose that R C g+ zs bounded set. Consider the functional 

/Cr(<p|t) 

wrt/z fl boundary condition r satisfying \\T\\oo,d"R < "/e- Then Kr has a unique maximizer 
{i^x)xGR and 

\\^\\oo,R < ||l'||oo,3"R- 

Proof of Lemma 9.5. Clearly we may assume £ [—"/2, ^/i], otherwise the configuration 
cannot maximise )Cr. Suppose first that is a maximizer and that t?x G ( — ^/2, ^/i) for 
all X e R. 

We may think of the equation satisfied by a stationary point of ICr in the form 

- E C-^i/ i^- - ^y] + = 0- (9-10) 

J/~.T 

where the conductances and potentials are given in terms of i? by 

sin(V(e,!?)) , „ 
C, = ^^A_i_i2 for alien R^0, 

= -m, cos(i?,)^^^^ for all x G Ql- 

Because we are assuming ^x G (-'^/2, '^/2) for all :c G R, Cg > < 0. 

Thinking of Cg, Vx as given and fixed, it is natural to define corresponding linear oper- 
ators on functions / on R U 3°R. 



• A = ^ E C^ytA - /y] for X G R,y G R U a°R 

'^X y^x 



■fx = for X G R 

We may rewrite (9.10) as 

{-L^ + M^)-fx = (9.11) 
subject to the boundary condition /| 90 r = t. The operator {—L^ + M^)^^ is positivity 
preserving, so if denotes the solution with boundary condition r' then 

<r < and||/l"l||oo,R < |!T||oo,aoR. 

Hence, if we can show that maximizers do not take on the values { — "/i, ^/i} the lemma 
will be proved as K, is uniformly convex in the region 

dx e [-^/6,^/6] Vx G R. 
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To show this let us in introduce a vector field on [— "^/2, "^/2] x R of the form 



x{e,x) 



-[0 + ^/6] if0<-^/ 
otherwise. 



and consider the flow dt6{t,x) := X{6{t, x), x). Then from any initial configuration 
9x{0) e [—"/2, ^/i] we find that /C(0(t) |t) is increasing, strictly so if there is a space time 
point {t, x) so that x) | > Our claims then follow immediately. □ 



We next claim that the maximizer i?j of K.Q+ {6\t) is very small for any x G so that 
either dist(x,9°fl+) is (much) larger than e^^ or dist(x,M n [ttAn U ^l^^J) small enough. 

Lemma 9.6 (Bulk Behavior) There exists eo G (0,1) so that ifO<£<£o the following 
holds. Let d optimize 

/Cg+(z?|T) with ||T;i-||oo,a»g+ < "^/e- 
Then for any x G g+ so that dist(x, B^g^) > Vs 

< \dx\ < exp — clog^'e . 

Lemma 9.7 (Boundary Behavior) There exists eo G (0, 1) so that ifO<£< £q the following 
holds. Let d maximize /Cg+ ['&\t) subject to a boundary condition r with 

||tx|U,3°0+ < ""/e and = Ofor all x G a°g+ n [9" Am U ^t/J- 

Then 

< I^^.tI ^ exp [ — clog^e] 
whenever dist(x,M+ n U ^t/J) < ^U- 

Together these Lemmas imply 

Corollary 9.8 There exists G (0,1) so that if < e < £q the following holds. Ifdist{x,M) < 
V5 then 

0< <exp [-clog^e] 

The proofs require the following technical input. Let Xt be a continuous time random 
walk on Z'^ with conductances (Ce),,^^: and let be the probability measure for the 
process. 

Theorem 9.9 (Gaussian Bounds [6]) Suppose there exist constants A,B > so that the 
conductances A < Cg < B/or all e. Then there exist constants Ci, Q, Cg, Cg so that for all t > 1 
and all x G Z'' 

C^^-\-c,\\^-y&,A{i+\\^-y\\2) < QC^Xt = y]< qr'^^e-^.^"'"-'"''^'^(i+ll-^-J/!l2) (9.12) 
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Proof of Lemma 9.6. The proof uses the Feynman-Kac formula and Azuma's inequality. 
Let J? be the unique maximizer of /Cg+ and set 

_ sin(V(e,i?)) 



V(e,t?) 



for all e n 0+ 7^ 0, 



sin(!?j) 

Vx = —mj cos(t?x) for all x G g+. 

By Lemma 9.5, there are constants A, B > so that A < l^^l/,,,^, Cg < B. In order to be 
able to apply the bounds of Theorem 9.9, we extend Ce = 1 for all e such that e H g+ = 0. 

Let Xf be a continuous time Markov chain with conductances Ce started from x G g^, 
Qx denoting expectation of the corresponding measure. Let 

R = dist(x,9°g+) AL/2 

By assumption R > ^1$. We define the stopping time Vx{^/k) as the first exit time of Xt 
from {y : ||x — y||oo < ^/k} and set Vq+ to be the first exit of Xt from 0+. Using the 
Feynman-Kac representation and the a priori bound on \ Vx\, 



-Jo' \Vx,\dt 



mx,df 



(9.13) 



Let Tt denote the natural filtration of Xf/^^, ^1^1/2). For any sequence of times 

So = < si < ■ • • < sj < dist(x,9°g+)^ = sj+i. 



let 



and note that 

v,{l/2) 





7+1 



s/+iA!;.,(l/2) 
SiAv,{l/2) 

J 



mx^dt - Qx 



s;+iAu,(l/2) 
s;Ai;,-(l/2) 



mx,dt\Tsi 



mx,dt > 3i + J^l{si < Ux-(l/2)}Qx,, / ' mx,dfl{s,+i - < Vx{l/2)} 

where the strong Markov property has been used to evaluate the second term and we 
note that the first sum is of martingale increments. We choose s,+i — s, = log^^°(L) for 
all /. 

Recall that Q'^ denotes the set of shifted cubes {Q(r + ^) : Q E Qi, rj E Af}. For each 
X E g^, we can find a cube Q of sidelength L so that 

QE Q[ and Ei{Q) = l, 

{y: ||x-y||oo <i^} C Q. 

Consider the increment 

I = Qz mx,dn{si+i - Si < Vx{l/2)} 
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By Lemma 9.4 and because C ^Zi{Q')=iQ' 

Because E(Q) = 1 the event -^jogso ^ occurs in the box Q. For e small enough, the Gauss- 
ian bounds in Theorem 9.9 imply 

I^e^[s/+i-Si]. 

as long as z G {y : ||x — y||oo < ^/i}- 
We have 



El{s,<i;v(V2)}Qx,. / mx,dtl{s,+i -s,- < t;;,.(V2)} 

!=0 



> X]l{s, < y,(V2),X,, G {y : |!x-y||eo < ^/i} e%+, - Si]. (9.14) 



! = 



Let M > 2 be fixed. Introducing the event F = {i;j(l/4) > ^^/m}, we thus obtain 



rVxn/2) 

-Jo "'Xfdt 



Here Ci, C2 are universal constants and the last term comes from applying the Gaussian 
estimate of Theorem 9.9 to the event f 
Using Lemma 9.4 once again, 

|V(e,gJ'?)l<elog^«e, 



so that 



240 



where we used log L ~ | log e| . 

Because the 3, are martingale increments, Azuma's Inequality then provides the bound 



< c/[.2log^«(.)]2 



Gathering the estimates together, we have proved that 

< tx < exp [ci (-^Vm + e^) L^] + exp [-C2M] . 
For e small enough, optimizing with respect to M proves the Lemma. 



(9.15) 

□ 



Proof of Lemma 9.7. We begin with the following observation. Using notation from the 
previous proof, suppose that x S 0+ and dist(x,M n [9' An U ^^l^^^J) ^ 
We may write 



\^x\ < Qx 



(9.16) 
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where Vg+ denotes the first exit time of 0+ and A < l^^l/m,- If we let 

then because we have assumed that the boundary condition Tx — for all x G 9° An U \0^]/ 



Let Br{x) denote the t'^ ball of radius r around x. Then if x G and 

dist(x,M n [9° An U '^t/^)) < 

If fx denote the first exit time from Bl/^(x), 



|l?rl < 



However a small variation on the argument in the previous proof gives 



and so the bound follows immediately. 



□ 



10. Proofs for §§ 5, 6 and 7 
10.1 Proofs for § 5. In this section we will prove Lemma 5.1 of § 5. 

Proof of Lemma 5.1. From our construction, the only statement which requires further 
justification is the first. Consider first the contributions from the interactions between 
neighboring boxes Qt^^{xo),Q£/^{xi) G Qf/^{T). For ANY edge e = {x,y) where x G 

Q%(^o), y e Q£/,(xi) for xq + xi, \V,a^(^^\ < e\ logep^. Thus 

< l-(Tx-(Ty<l. 

e=(x,y): 

As a consequence, our choice of I implies 

0< E MQ) + ns(Ty^'^)<Ce"no^''\T\+ E E(Q) + 7^Q(^"(r))_ 

where we used the fact that E(Q) > Eo(Q) on the RHS. 

Fix Q such that H(y^(Q) = 1. In this case, a combination of Proposition 8.1 and (8.5) 

implies that \£Q{g^) is (essentially) the maximum value —T-Lq can be on Sq. Consider 
the free boundary condition maximum 

E(Q)=max--HQ((7) 
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We have 

" . 

I 



+ 



-nQM-^£Qi§Q) 



(10.1) 



n III 

N 



with (TQ^y was defined at (5.1). Examining the term ErrovQ in Lemma 8.2 and recalling 
the conditions which entail E/ , = 1 

'2 

II + III < el 

with plenty of room to spare for e small. Regarding Term I, using that \t9\ < 1 and 
Ze/^{Q) = 1, if we only display dominant terms involving r, we have 

I < II^^IILII Vrlli^Q + II VgglloollV^glblll - T||2,Q + el 
Let W = {:c G Q : Tx 7^ 1}. Clearly 



I-t^IIIq.. <\W\<£-'H^. 



and 



Using the bounds H^'qUoo < 2e£2|logep° and ||V^q^, ||oo < e|loge|'^'^. 



1 30 



I < elt. 



Thus we have 



Also, 



^ [e(Q) + UQ^a'^H 1{S.,^(Q) = 1} < el|r|. 



E [E(Q) + 7^Q(^r^"(r))jl{H,,^(Q)=0} 

< E(Q)l{S,,^(Q) = 0,fo^(Q) < e^llogel} 



X: E(Q)l{E,,^(Q)=0,fo^(Q)>e2|loge|} 



where we used that 'Hq(c7^(^') = on these boxes. Using (3.14) and (8.1) and the fact 
that the density of boxes with Ht^ = is small to bound term I and using (3.11), (3.14) 
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and (8.1) to bound term II 

I + II < e2|ioge|-40|r|. 

□ 



10.2 Proofs for § 6. In this section we prove Lemmas 6.4, 6.6 and 6.7. For notational 
convenience, we will use the functions Fa/ F^' ^ defined after (3.10) in our estimates 
with A = L-^log^L. 

Proof of Lemma 6.4. We restrict the discussion to Using (8.3) along with the fact that 
cr^ = on 9°i^+ to cancel boundary contributions, we have 



where ^' = QJ^ U 3°^+. By definition, the original spin configuration a has '^^^^ =1 on 

e:(r) so 

Eoj\a)<e^\\o^e\\^^\. 

However if if F is clean. 



^ |loge 



-55 



by condition (3.10). Thus Eoj,{(j) < e^l loge|-50 
Next, we have by the Cauchy-Schwarz inequality 



2 



QgSl QgSl 
and by (3.13) 



QeSL 



Therefore, by Lemma 9.4, we have, for e small enough, 

^i^'(4'?)<i' E Fl +e-'^'°^'^'n (10.3) 

^ V ' 

I 

QgQl 

We bound the first term on the RHS of (10.3) via 

I<e2|ioge||^+|+L'* J2 FA^(Q)l{fA^(Q)>e'|loge|}, 



QeQi 



<e'\^oge\-^\T\ 
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To obtain the second bound, we used the fact that | . 
A and (3.11) to bound B. 

A similar estimate on \\g^^ Wi using (3.12) in place of (3.11) gives 

Combining the estimates yields the claim in the Lemma. 
Proof of Lemma 6.6. Note first of all that by construction. 



< \D\ < |loge| ^^|r| toboimd 



□ 



where a is the original spin configuration from which cr" is produced. Therefore, Lemma 8.2 
implies 



-nA^{o^\ei) > -HaM^i) - CiError^^ (a) + Error! 



Cie I log e 



-25 I 



Because Hl(Q) = 1 for {Q G Q'^^ : Q n 'r+ 7^ 0}, inequalities (3.8) and (9.7) allow us 
to exchange terms involving g^i^ for sums involving {gQ^)Qnv+^i!i up to errors of order 

g-ciog*e_ Thus 

Error^J',{a) + Error^^,{g^) < Error^Q^" {cr) + Error^f {g'f) + e^\\o^e\-^^\T\ 

Q^Q'l 

Recall that Sqicr) < \ loge||Q| if Q n "tf^ ^ because Q is a cube at the outer 
boundary of a contour. Then because ^l{Q) = 1/ inequalities (3.4), (3.5), and (3.6) imply 

£rror^'°((7) + Error'f {g'^'' ) < e^^\Q\ (10.5) 

for £ small enough and the claimed lower bound follows. □ 



Proof of Lemma 6.7. To begin with 



< 



V 

I 



+ 



(10.6) 



To estimate these terms we need some bounds on the interpolation. First, let 

W = {x:3y~x,c7-yV^^} 
denote the support of the modification. Clearly 

|W| < i-'^\Y\ and 
For Term I, we have, using corollary 9.8, 



(10.7) 

(10.8) 
(10.9) 

(10.10) 
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For Term II, write eu = [—A + A] • g^^on ^. Summation-by-parts gives (in particular 
using that since W C is strict to cancel boundary terms) 



II < 



eCW 



+2A E l^'Sl 
xew 



III 

We then have, using corollary 9.8, 

iii<£:w((73) + ^^(g.^'^)+ E ^'f [V(.,)T]2 + e-ii°g^iVl- 

{xy)cW 

We now derive bounds on the respective RHS's. By corollary 9.8 and the definitions 
of W, (P, 

£w{(J^) < Sw{g^f)+e-''°s'e\Y\. (10.11) 
Since ^ C U2j^(q)=iQ, Lemma 9.4 implies 

aE l^^'5l<e''^|r|, (10.12) 

Sw{g^f)<e'*\T\. (10.13) 

because llgg'^llco < e^/Llog^^eand ||V^q'^||oo < e|logepo if Si.(Q) = 1 and also using 
(10.8). Similarly, using (10.8) and (10.9) 

{xy)cW 

and the claim is proved by collecting these estimates. □ 

10.3 Proofs for § 7. Recall that we defined = {Q(r) : Q(r - (^/i, ^/i, ^li) e Q{}. 
Proof of Lemma 7.1. Let (7 G X(r) and set 

T^eiT) = {QeQr.Qn 5(1) ^ 0}, 7^^(^) = {QeQ',:Qn <5(r) ^ 0}. 

We first argue that enough of the energetic defect of a is captured by restricting attention 
to its behavior to boxes either in 7?.^(r) or in 7?.^ (F). 
Let 

5^ (F) = {Q C F such that Q e Ql and Y;,^ (Q) = because ^i^^ = for some z}, 
5Ave(r) = {Q C F such that Q E Ql and Y;,(,(Q) = but ipi^^ = 1 for all relevant z}. 

If Q G 5^ (F) there exists a cube of side-length £, Qi, so that dist(Q^, Q) < 2L + 5£ and 

%,((7) > epilog e|£3. 

Then there must be Q' in either 7?.£(F) or 7?.'^(F) so that 

%(t7) > <^Vi6|loge|£3. (10.14) 



42 N. CRAWFORD 

Note that these considerations include cubes which overlap in the case of 7?.^(r). Let 

Ai = {Q G 7^^(^) such that (10.14) holds} (10.15) 

A2 = {Q G n'f{Y) such that (10.14) holds} (10.16) 

On the other hand if Q G SAvei^) there exits so that dist(Q/v Ql(^o)) < 2L + 5i 
and 

k(Qf)-ei| < 1-?. 

Then there must be Q' in 7?.c (F) 

\(r{Q')-ei\<l-^ (10.17) 

Let 

A3 = {Q G Te^r) such that (10.17) holds}. 

Because (7- G X(r), 

max|A,-| >\lose\-^^Nf^\ 

i 

Now we compare the internal energy of a with that of S^. Recall that ext denotes the 
boundary condition which is set to e\ on 9°R H and is free otherwise. Then reflection 
invariance of the Hamiltonian for components with free boundary conditions implies 

-'^AM\J(r)(t^^|ext) = -'H^^\j(r),(t7*|ext). 

Thus, by Lemma 6.1 

-7^AjS±|ei) +7^Aj^ki) > -7^j(r)(c^'("^|ext) +7^,-(r)(c7'^|ext) -Ce^l logel^ 

^ ^/ ' 

I 

(10.18) 

Next recall that 

Eo(Q) = max -HQ{a\ex\.). 

By Lemma 5.1, 

-7^,(r)(c.^"(i')|ext) > £ Eo(Q) - e^l loge|-40|r| 

Note that if we have a finite region R = Ri U R2 with Ri n R2 = then 

Eo(R) < Eo(Ri)+Eo(R2). 
Note also that if Q G A;, Q n A^j may be covered by boxes in TZc/ (F). Therefore, for 

2 e {1,2,3}, 

I> ^ [Eo(QnAN)+-HQnAMMext)]-e2|loge|-25|r| 

QeA; 
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The superscript C was dropped because the original spin configuration a is unmodified 
on such boxes as they must intersect sp(r) . It remains to show that 

maxll, > e2|ioge|-24|r|. (10.19) 

i 

We consider only A2 in detail. The remaining cases are a bit simpler to handle as we 
do not have to deal with cubes overlapping A^. 

Since we can always choose (Ty = e\, Eo(Q n A^f) > 0. If Q G A2 and ^^(Q) = \, (8.1) 
of Proposition 8.1 implies 

--HQnA.Mext)<-e2|loge|£3 
because we may view — 'HgnAwlc^lext) as —'HQ{a') for the configuration 

^, _ /c7j/if y e Qn An 
ei ify G QnA^. 

Thus 

^ [Eo{Q)+nQ{a\ext)]>e^\loge\i'\A2\. 

QeAr.EeiQ)^! 

On the other hand, using (8.1) once again, 



X: [Eo(Q)+?^Q((7|ext)]>- X] SQig^) + e\^m^' 

QgA2:SKQ)=0 QeA2:S((Q)=0 



III 

Using conditions (3.11) and (3.14) and that fact that S{Y) is good at the scale i, 





ni < 


loge 


-50 p 


It follows that 








112 > 


\oge\f 




e^\\oge\-''\] 


The same sort of argument gives 








ni>^^| 


loge|£3 




e^\loge\-^°\ 



(10.20) 

II3 > ?W|A3| -e'lloger^^llrl (10.21) 

where we have used (8.5) instead of (8.1) to estimate the contribution from A3. (10.19) 
and hence Lemma 7.1 follow. □ 



11. Estimates on the Randomness 

In this section we derive elementary probabilistic estimates which underly the rest of 
the paper. Let I G N and restrict a to Q/ C Z"^. With an eye toward future work we 
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record bounds for all cf > 2 though we only use the d = 3 case here. 

(-Ag,+A)-i.a, 











Ql ~ 






if 






Ql ~ 


£ 



(-A~ +A)-i-[«-*Q,]. 

Here AG [0/1). From now on Q/ is fixed and we drop this subscript from our notation. 
It is important for us to have fairly precise probabilistic estimates on the quantities 

VU oo/ oo, 2/ VLr 2/ 



and similarly for G"^'^, to have estimates on ocq^ \/|Q/| and probabilistic bounds on fluc- 
tuations of the low momentum modes. In d = 2 we record somewhat more refined 
information, in particular bounds on the density of points with atypical fluctuations. 
For general finite regions R C Z'^, we introduced potentials of the form 

xee 

in § 8 and used its typical behavior heavily in § 6. Here we specialize R = Q/. For 
notational convenience set 

R^{x) = A-3 Adist(x,a°Q;) 
Given A > 0, consider the event 

Ar = Ar{A, Q,) = J a; : r'^ ^ nty > Ae^log^^-'' R^{x) for all x e Qi s.t. distoo(x,a°Q/) 

[ ||y--i;|U<r 

where S2^d is 1 if = 2 and otherwise. 
Let 



_[ \\k\\^ 



The following lemma summarizes the bounds we need. 

Lemma 11.1 Let d > 2 and 5 G (0,1) be fixed. For any A G [0,1) we have the following 
probabilistic estimates for either choice of boundary conditions: 



(1) Let M G (1, oo) be fixed. For any x G Q/ and any edge e 

p(|G,^|>MG2)<e-^^' 
p(|V.G'^|>MGv)<e-^^'. 
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(2) 



w{\\GY2>m(;\n < 



/-™^ for d = 4, 

log log /g-"'" /or d G {2,3}. 



(3) For all d >2, there is > so that 



F{\\VG^\\l<c,gll')<e 



(4) 



,-cml'' 



, for d > 3, 



log log /e""" ford = 2. 



F(^\\VG'\\l>mgll')<< 
(5) There is > so that ifO<\< /"^ 



We will say a bit in the way of proof about these bounds in a bit, but let us first address 
Proposition 3.3 and Lemma 3.4. 



Proof of Proposition 3.3. We show 



P 



U-'^A 'l{FV>.2|iog£|}J^ >be Ny j ^e 

only. All remaining bounds are proved in a similar or simpler way using Lemma 11.1. 
For each;/ G {-32, . . . ,32}3 let 

T, = {Q{r + fj):QcY,QeQL,} 

for some choice oitj e if {-32, . . .,32}*^. Let 

m{Q) = min{m G N : HV^qHj < me^| loge|Lo^} 

To obtain the estimate, it is enough to bound probabilities for the finite collection of 
events indexed by rj 

^ m(Q)l{m(Q) > 2} > e'^^N^" 

For any subset ^ C 7-",^, any choice {n(Q) : Q G A,n{Q) > 2} and for e small enough, 

P(mQ = n(Q) VQ E A) < d'^le-'^Oi^^"^^^^''^'^^^ 

This implies, after some standard computations to take into account the entropy of the 
family of subcubes which contribute to the sum, that 

P I ^ m{Q)l{m{Q) > 2} > MN^° ] < c^y° e-'^\^°&'\^o^Y° (ll.l) 
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for any M > 0. Taking M = e'"^ the exponent on the RHS is still at least e and the 

claim follows. 

□ 

Proof of Lemma 3.4. We claim first of all that 

T{A{x)) < Ce-'\^°S^\' 

and 

|P(A(x)nA(y))-P(A(x))P(A(y))| < Ce-'^°s'^distdx,y) (ij^) 

Here distL(x,i/) denotes the minimal number of blocks in an L-measurable block path 
from Q(x) to Q{y). 

Consider the collection of bounded L-measurable connected subsets of Z'^ containing 
X, {y}Y3j- For each such Y, we can apply the probabilistic estimates of Proposition 3.3. 
The number of L-measurable connected sets Y containing x with Ny = r is well known 
to have the asymptotic Aq for some fixed, dimension dependent constant uq. Thus 

P(Q(x) is in some Y which is not clean) < C ^(2flo)''e"'^°s'e'- if = 3 

)->i 

Modifying this estimate slightly via the discrete isoperimetric inequality, 

P(A(x) is in c(y) for some y which is dirty) < C ^ r'^/('^-^'(2flo)'V'°*^'^'' ifd = 3. 

(11.3) 

The first claim follows. 

Next we prove the correlation bound. Using the fact that the events {Y,- is dirty} are 
independent if (52l(c(Yi)) n^2L(c(Y2)) = we have 

|P(A(xi),A(x2)) -P(A(xi))P(A(x2))| < 

F(Q(^i), Qixi) are in c(Yi),c(Y2) for some dirty Yi, Y2 with ^2l(c(Yi)) n ^2l(c(Y2)) = 0). 

(11.4) 

An estimate similar to (11.3) then gives 

|P(A(xi),A(x2)) -P(A(xi))P(A(x2))| < c|Qi.|e"'^°s'<^dist,(^i,x2)_ (^j ^) 

The first bound implies 

E[|Dm|] <|AN|e-^|i°S'^l' 
The correlation bound implies that 

Var[|DM|] <L3|ANk-^l'°s^l' 

By taking N = 2^, applying Chebyshev's inequality to estimate deviations of \T)n\ and 
then the Borel-Cantelli lemma along this subsequence, we have that, for almost every 
a; G Q, there is No(a;) G N so that 

|An| ~ 
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□ 

Proof of Lemma 11.1 (1),(2), (3), and (4). These bounds are elementary computations and 
we only sketch the basic argument. Statement (1) simply relies on the fact that G^, V^G'^ 
are Gaussian variables with mean zero and variance bounded above by a constant mul- 
tiple of g2, Gv respectively. The idea of (2), (3) and (4) is to expand the field G'^ in terms 
of either Dirichlet or Neumann Laplacian eigenfunctions depending on the boundary 
conditions. By (a relative of) Parseval's identity, we may pass to momentum space and 
express all quantities of interest as weighted sums of squares of i.i.d standard Gauss- 
ian variables (i.e. weighted by the eigenvalues of [—A + A]^'^, — A[— A + A]^^). Here is 
where, for convenience, we use the Gaussian assumption: Fourier transforms of i.i.d. 
Gaussians are i.i.d. Gaussians. 

These latter sums are estimated by first separating summands according to the mo- 
mentum space annuli 

As := {k G ...l}^: 2-(^+i) < ||ic||2 < 2-^} 

This is useful because eigenvalues of A corresponding to these momenta are the same 
up to a multiplicative constant independent of s and we can treat the contribution from 
each annulus as a constant multiple of an i.i.d. sum of squares of Gaussians indexed 
by wave vectors in the annuli. If the cardinality of As is big enough, the corresponding 
sum is highly concentrated around its mean while if the cardinality of As is small the 
corresponding sum of squares must have a reasonably large fluctuation to contribute to 
the overall summation. The extent to which these sets really contribute is reflected in 
the various cases stated in the Lemma. □ 

Proof of Lemma 11.1 (5). Let us define 

n,- = {Q C Q/ : Q is a cube of sidelength r}, (11.6) 

F = n,>iog«o, DQ^n.. : liVg^'^lli > Cde\l^,/}, (11.7) 

fi = {a; : || V.g^f ||oo,Q,uaoQ, < elog'^ I}. (11.8) 
Then by Lemma 11.1 (1),(3), 

P(ruF0<exp(-clog55Z). 

We prove that f n Fi C Ar{Aij) for appropriate choices of r,Aci > 0. On F n Fi, let 
Q e ^r>ios'° i^'-- X E Q,we may express the field g'^^^ via 

„A,D A,D I ^(1) 

SQi,x = 8Q,.r 

where g^^^ is satisfies the Laplace equation —Ag^^^ = on Q subject to the boundary 
condition g[^^ = g^^^ for x E d°Q. 
Notice that 

E mg'£+g^^^)f = £Qu3^Q{g'£)+ E [Ve^i'^]' 

enQ^0 enQ/0 
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since the cross term vanishes. This is because g^^^ is harmonic in Q and gQ^ vanishes on 
a°Q. Hence 



90 



because we restricted attention to f n fi. Thus f n f i C Ar{cd, Qi) whenever log / < 
r < '/4. □ 
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